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Abstract 

A sample of 114 ± 11 B® — > J/ipK~n + signal events obtained with 0.37 fb _1 of pp 
collisions at yfs = 7 TeV collected by the LHCb experiment is used to measure the 
branching fraction and polarization amplitudes of the B® — > J/tpK*° decay, with 
K*° — > K~n + . The K~ir + mass spectrum of the candidates in the B® peak is dom- 
inated by the K*° contribution. Subtracting the non-resonant K~ir + component, 
the branching fraction of B® — > J/ipK*° is (4.4 j^Q \ ± 0.8) x 10~ 5 , where the first un- 
certainty is statistical and the second is systematic. A fit to the angular distribution 
of the decay products yields the K*° polarization fractions fi = 0.50 ± 0.08 ± 0.02 
and /|| = 0.19±g;Jg ± 0.02. 
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Interpretations of measurements of time-dependent CP violation in B® — > J/ip <fi and 
B° s — > Jjijj / (980) decays have thus far assumed the dominance of the colour-suppressed 
tree-level process. However, there are contributions from higher order (penguin) processes 
(see Fig. [T]) that cannot be calculated reliably in QCD and could be large enough to 
affect the measured asymmetries. It has been suggested that the penguin effects can be 
determined by means of an analysis of the angular distribution of B® — > J /ipK* (892)°, 
where the penguin diagram is not suppressed relative to the tree-level one, and SU(3) 
flavour symmetry arguments can be used to determine the hadronic parameters entering 
the B° s —7- J/ip(j) observables [l]. 

In this paper the i^*(892)° meson will be written as K*°, while for other K* resonances 
the mass will be given in parentheses. Furthermore, mention of any specific mode implies 
the use of the charge conjugated mode as well, and K~n + pairs will be simply written 
as Kit. The decay B Q S — > J/ipK*° has already been observed by the CDF experiment J2], 
which reported B{B° S — > J/ipK*°) = (8.3 ± 3.8) x 10~ 5 . Under the assumption that 
the light quark (s,d) is a spectator of the b quark decay, the branching fraction can be 
approximated as 

B(B° S ->■ J/^K*°) ~ x B(B° ->■ J/^K*°) = (6.5 ± 1.0) x 1(T 5 , (1) 

I ' cs I 

with \V cd \ = 0.230 ± 0.011, \V CS \ = 1.023 ± 0.036 (3), and B(B° J/i)K* Q ) = (1.29 ± 
0.05 ± 0.13) x 10 -3 |4j. The measurement in Ref. |4j, where the Kir S-wave contribution 
is subtracted, is used instead of the PDG average. 

In this paper, 0.37 ftT 1 of data taken in 2011 are used to determine B{B° S — > J/i/jK*°), 
to study the angular properties of the decay products of the B° s meson, and to measure 
the resonant contributions to the Kir spectrum in the region of the K*° meson. The 
measurement of the branching fraction uses the decay B° — > J/ipK*° as a normalization 
mode. 

The LHCb detector |5] is a single-arm forward spectrometer covering the pseudo- 
rapidity range 2 < i] < 5. The detector includes a high precision tracking system con- 
sisting of a silicon-strip vertex detector located around the interaction point, a large-area 
silicon-strip detector located upstream of a dipole magnet with a bending power of about 
4Tm, and three stations of silicon-strip detectors and straw drift tubes placed down- 
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Figure 1: Tree and penguin decay topologies contributing to the decays B® — > J/ipK*° and 
Bg — > J/tp(j). The dashed line indicates a colour singlet exchange. 
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stream. The combined tracking system has a momentum resolution Ap/p that varies 
from 0.4% at 5GeV/c to 0.6% at lOOGeV/c. Two ring-imaging Cherenkov detectors 
(RICH) are used to determine the identity of charged particles. The separation of pions 
and kaons is such that, for efficiencies of ~ 75% the rejection power is above 99%. Pho- 
ton, electron and hadron candidates are identified by a calorimeter system consisting of 
scintillating-pad and pre-shower detectors, an electromagnetic calorimeter and a hadronic 
calorimeter. Muons are identified by alternating layers of iron and multiwire proportional 
chambers. 

The trigger consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage called High Level Trigger (HLT) that 
applies a full event reconstruction. Events with muon final states are triggered using two 
hardware trigger decisions: the single-muon decision (one muon candidate with transverse 
momentum p^ > 1.5GeV/c), and the di-muon decision (two muon candidates with pr,i 
and pt,2 such that y/pT,iPT,2 > 1.3GeV/c). All tracks in the HLT are required to have 
a pt > 0.5GeV/c. The single muon trigger decision in the HLT selects events with at 
least one muon track with an impact parameter IP > 0.1 mm with respect to the primary 
vertex and > 1.0 GeV/c. The di-muon trigger decision, designed to select J/ip mesons, 
also requires a di-muon mass (M^) 2970 < M m < 3210 MeV/c 2 . 

Simulated events are used to compute detection efficiencies and angular acceptances. 
For this purpose, pp collisions are generated using Pythia 6.4 [6] with a specific LHCb 
configuration |7|. Decays of hadronic particles are described by EvtGen [8] in which final 
state radiation is generated using Photos j9]. The interaction of the generated particles 
with the detector and its response are implemented using the Geant4 toolkit 10 as 



described in Ref. 11 



The selection of B^ — > J/ipK*° decays first requires the reconstruction of a J/ip — > 
candidate. The J/ip vertex is required to be separated from any primary vertex 
(PV) by a distance-of-flight significance greater than 13. Subsequently, the muons from 
the J/ip decay are combined with the K and ir candidates to form a good vertex, where the 
di-muon mass is constrained to the J/ip mass. A p^ > 0.5 GeV/c is required for each of the 
four daughter tracks. Positive muon identification is required for the two tracks of the J/i/j 
decay, and the kaons and pions are selected using the different hadron probabilities based 
on combined information given by the RICH detectors. The candidate momentum 
is required to be compatible with the flight direction as given by the vector connecting 
the PV with the candidate vertex. An explicit veto to remove B + — > J/tpK + events is 
applied, as they otherwise would pollute the upper sideband of the B9, mass spectrum. 

Following this initial selection, several geometrical variables are combined into a single 
discriminant geometrical likelihood variable (GL). This multivariate method is described 
in Refs. |l2|l3| . The geometrical variables chosen to build the GL are: the B^ candidate 
minimum impact parameter with respect to any PV in the event, the decay time of the 

v candidate, the minimum impact parameter x 2 of the four daughter tracks with respect 
to all PV in the event (defined as the difference between the \ 2 °f the PV built with and 
without the considered track), the distance of closest approach between the J/ip and K*° 



2 



trajectories reconstructed from their decay products, and the px of the B9, candidate. 
The GL was tuned using simulated B° — > J/tpK*° signal passing the selection criteria, 
and background from data in the B^ mass sidebands with a value for the kaon particle 
identification variable in a range which does not overlap with the one used to select the 
data sample for the final analysis. 

The K7i mass spectrum in the B° — > J /if) Kit channel is dominated by the K*° reso- 
nance but contains a non-negligible S-wave contribution, originating from _K"q(1430)° and 



non-resonant Kti pairs 14 . To determine B(B° S — > J/ipK*°) it is therefore important to 
measure the S-wave magnitude in both B9, — > J/tpKrc channels. The Kit spectrum is 
analyzed in terms of a non-resonant S-wave and several Kir resonances parameterized us- 
ing relativistic Breit-Wigner distributions with mass-dependent widths, following closely 
[li] . The considered waves are: a non-resonant S-wave amplitude interfering with the 
i^o(1430)° resonance, K*° for the P-wave and fQ(1430) for the D-wave. F-wave and 
G-wave components are found to be negligible in the B° fit. In bins of the Kir mass, a 
fit is made to the B9s candidate mass distribution to determine the yield. As shown in 
Fig. [2J a fit is then made to the B° and B® yields as a function of the Kir mass without 
any efficiency correction. The S and P-wave components dominate in the ±40 MeV/c 2 
window around the K*° mass, where the K*° contribution is above 90%. A more exact 
determination of this contribution using this method would require Kir mass-dependent 
angular acceptance corrections. For the branching fraction calculation, the fraction of K*° 
candidates is determined from a different full angular and mass fit, which is described next. 

The angular and mass analysis is based on an unbinned maximum likelihood fit which 
handles simultaneously the mass (Mj^k-k) and the angular parameters of the B9 \ decays 
and the background. Each of these three components is is modelled as a product of 
probability density functions (PDF), V(Mj/^ Kn ,ip,9,{p) = V(Mj/^ K7T )V(ip,9,{p), with i[> 
the angle between the kaon momentum in the rest frame of the K*° and the direction 
of motion of the K*° in the rest frame of the B. The polar and azimuthal angles (6, (p) 
describe the direction of the /i + in the coordinate system defined in the J/ip rest frame, 
where the x axis is the direction of motion of the B9 \ meson, the z axis is normal to the 
plane formed by the x axis and the kaon momentum, and the y axis is chosen so that the 
y component of the kaon momentum is positive. 

The function describing the mass distribution of both B9, signal peaks is the sum of 



two Crystal Ball (CB) functions 16 , which are a combination of a Gaussian and a power 



law function to describe the radiative tail at low masses, 

V{M J/i)KlT ) = f CB(M J/i!K7r , fi B , 01, ai) + (1 - /) CB(M J/ ^^ 7r , fi B , a 2 , a 2 ). (2) 

The starting point of the radiative tail is governed by a transition point parameter 0(1,2) • 
The mean and width of the Gaussian component are [is and <J(\$)- The values of the /, 
Oi, o" 2 , a.\ and a 2 parameters are constrained to be the same for the B® and B° peaks. 
The difference in the means between the B® and the B° distributions, ([Ib° — Hb°), is 



fixed to the value taken from Ref. 17 . The mass PDF of the background is described by 
an exponential function. 
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Figure 2: Fit to the Ktt mass spectrum for (a) B° — > J/ipKn events, and (b) — > J/ipKir 
events. The £?9n — > J/ipKir yields in each bin of Ktt mass are determined from a fit to the 
J/ip Kit mass spectrum. The pink dashed-dotted line represents the K*°, the red short-dashed 
line is the S-wave and the black dotted line is the -^2(1430). The black solid line is their sum. 



Assuming that direct CP violation and the B9 S n — B^ production asymmetry are 
insignificant, the differential decay rate is [1,15 



d 3 r 

-j^y oc 2 1 Aq 1 2 cos 2 tp(l — sin 2 9 cos 2 <£>) 

+ |A||| 2 sin 2 ^(l - sin 2 6» sin 2 y?) 
+ |A ± | 2 sin 2 ?/>sin 2 6) 

+ |A | I cos(<5|| — <5 ) sin 2?/^ sin 2 6 1 sin (3) 
V 2 

+ l\A s \ 2 [1 -sin 2 #cosV] 
4^/3 

+ 1 A 1 1 As I cos(5s — S ) cosip [l — sin 2 #cos 2 ip\ 

3 

\/6 

+ — -|A|| I \A$\ cos(<5|| — 5s) sin-0 sin 2 6 sin 2<p, 

where v4 , Ay and A± are the decay amplitudes corresponding to longitunaliy and trans- 
versely polarized vector mesons. As = \As\e lSs is the Kn S-wave amplitude and (8\\ — So) 
the relative phase between the longitudinal and parallel amplitudes. The convention 
Sq = is used hereafter. The Q differential is dQ = dcos^dcos#dv?. The polarization 
fractions are normalized according to 

f 11 = I^o,ii,-lI 2 U) 

JL ' lU \A \ 2 + |A||| 2 + \A ± \ 2 ' y } 

which satisfy ft + f\\ + f± = 1. 
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The parameters fi, /n and 8\\ describing the P-wave are left floating in the fit. The 
\A S \ amplitude and the 5s phase depend on M Kn , but this dependence is ignored in 
the fit, which is performed in a K7r mass window of ±40MeV/c 2 , and they are just 
treated also as floating parameters. A systematic uncertainty is later associated with 
this assumption. The angular distribution of observed events is parameterized as a 
product of the expression in Eq. [3] and a detector acceptance function, Acc(f2), which 
describes the efficiency to trigger, reconstruct and select the events. Simulation stud- 
ies have shown almost no correlation between the three one-dimensional angular accep- 
tances Acc^ip), Accg(9) and Acc v (<p). Therefore the global acceptance factorizes as 
Acc(f2) = Ace,/, (ip) Acc,j(6>) Acc v (y2), where Acc^ip) is parameterized as a fifth degree 
polynomial, Acce(6 l ) as a second degree polynomial and Acc</,(</>) as a sinusoidal function. 
A systematic uncertainty due to this factorization hypothesis is later evaluated. The an- 
gular distribution for the background component is determined using the upper sideband 
of the B® mass spectrum, defined as the interval [5417,5779] MeV/c 2 . 

Figure [3] shows the projection of the fit in the Mj/^ Kn mass axis, together with the 
projections in the angular variables in a window of ±25 MeV/c 2 around the B® mass. The 
number of candidates corresponding to B° and B® decays is found to be 13,365 ± 116 and 
114 ± 11, respectively. 



Table 1: Summary of the measured — > J/ipK*° angular properties and their statistical and 
systematic uncertainties. 



Parameter name 


l^sl 2 


h 


h 


Value and statistical error 


0.07^; 5 7 


0.50 ±0.08 


0-19^ 


Systematic uncertainties 


Angular acceptance 


0.044 


0.011 


0.016 


Background angular model 


0.038 


0.017 


0.013 


Assumption 5s(M Kn ) = constant 


0.026 


0.005 


0.002 


B° contamination 


0.036 


0.004 


0.007 


Fit bias 






0.005 


Total systematic error 


0.073 


0.021 


0.022 



Tables 

ters, toget 



and 2 summarize the measurements of the B9 g ^ — > J/ipK*° angular parame- 
ter with their statistical and systematic uncertainties. The correlation coeffi- 



cient given by the fit between fi and f\\ is p = —0.44 for B® decays. The results for the 
B° — > J/ipK*° decay are in good agreement with previous measurements [4| [l5|[l8|[T9] . 
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Figure 3: Projections of the fit in Mj^ Kw and in the angular variables for the mass range 
indicated by the two dashed vertical lines in the mass plot. The red dashed, pink long-dashed, 
and blue dotted lines represent the fitted contributions from B° — > J/ipK*°, B® — > J/ipK*° and 
background. The black solid line is their sum. 



Based on this agreement, the systematic uncertainties caused by the modelling of the 
angular acceptance were evaluated by summing in quadrature the statistical error on 



Table 2: Angular parameters of B° — > J/ipK*° needed to compute B(B® — > J/ipK*°). The 
systematic uncertainties from background modelling and the mass PDF are found to be negligible 
in this case. 



Parameter name 


l^sl 2 


h 




Value and statistical error 


0.037 ±0.010 


0.569 ± 0.007 


0.240 ±0.009 


Systematic uncertainties 


Angular acceptance 


0.044 


0.011 


0.016 


Assumption 5s(Mktt) = constant 


0.026 


0.005 


0.002 


Total systematic error 


0.051 


0.012 


0.016 
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the measured B° — > J/tpK*° parameters with the uncertainties on the world averages 
(/l = 0.570 ± 0.008 and f± = 0.219 ± 0.010) pi. The angular analysis was repeated with 
two additional acceptance descriptions, one which uses a three-dimensional histogram to 
describe the efficiency avoiding any factorization hypothesis, and another one based on a 



method of normalization weights described in Ref. 20 . A very good agreement was found 
in the values of the polarization fractions computed with all the three methods. For the 
parameter |^4s| 2 > uncertainties caused by the finite size of the simulation sample used for 
the acceptance description, as well as from the studies with several acceptance models, 
are included. The systematic uncertainty caused by the choice of the angular PDF for 
the background is shown for the B® — > J/ipK*° decay but it was found to be negligible 
for B° J/ipK*°. 

Also included in Tables [T] and [2] is the uncertainty from the assumption of a constant 5$ 
as a function of M^. This assumption can be relaxed by adding an extra free parameter to 
the angular PDF. This addition makes the fit unstable for the small size of the B° s sample, 
but can be used in the control channel B° — > J/ipK*°. The differences found in the B° 
parameters with the two alternate parameterizations are used as systematic uncertainties. 
The parameters 5\\ fit to cos(<J||) = -0.960i^ for the B° and to cos(5||) = -0.93 ± 0.31 
(where the error corresponds to the positive one, being symmetrized) for the B° s . These 
parameters could in principle affect the efficiency corrections, but it was found that the 
effect of different values of 8\\ on the overall efficiency is negligible. A simulation study of 
the fit pulls has shown that the errors on /t and h of the B° s decays are overestimated 
by a small amount (~ 10%) since they do not follow exactly a Gaussian distribution, 
therefore the decision was taken to quote an uncertainty which corresponds to an interval 
containing 68% of the generated experiments, rather than giving an error corresponding 
to a log-likelihood interval of 0.5. A slight bias observed in the pulls of f» in B° s decays was 
accounted for by adding a systematic uncertainty corresponding to 6% of the statistical 
error. 



The ratio of the two branching fractions is obtained from 

B(B° a -» J/jjK*°) _ f d eg \ B o /go N B o 



(5) 



where fd (/«) is the probability of the b quark to hadronize to B° (B®) mesons, e^J/eg is 

the efficiency ratio, \ B a/\ B o is the ratio of angular corrections, f^l /f^lo is the ratio of 
K*° fractions and N B o/N B o is the ratio of signal yields. The value of fd/ fs has been taken 



from Ref. 



21 



The efficiencies in the ratio e^l/e^l are computed using simulation and 
receive two contributions: the efficiency of the offline reconstruction (including geometrical 
acceptance) and selection cuts, and the trigger efficiency on events that satisfy the analysis 
offline selection criteria. The systematic uncertainty in the efficiency ratio is negligible 
due to the similarity of the final states. Effects due to possible differences in the decay 
time acceptance between data and simulation were found to affect the efficiency ratio by 
less than 1 per mille. On the other hand, since the efficiency depends on the angular 
distribution of the decay products, correction factors A^o and \ B o are applied to account 



7 



Table 3: Parameter values and errors for g(^oZ^j/^*o] ■ 



Parameter 


Name 


Value 


Hadronization fractions 


fd/ fs 


3.75 


± 


0.29 


Efficiency ratio 


^.tot / _tot 

t B o/fc B o 


0.97 


± 


0.01 


Angular corrections 


Abo/Abo 


1.01 


± 


0.04 


Ratio of K*° fractions 


As) i Ad) 
J K*°l J K*° 


1.09 


± 


0.08 


B signal yields 




(8.5183 ± 


0. 


8) x 10~ 3 



for the difference between the angular amplitudes used in simulation and those measured 
in the data. The observed numbers of B° and B® decays, denoted by Nbo and Nbo, 
correspond to the number of B® — > J/ipKn and B° — > J/ipKir decays with a Kit mass 
in a ±40 MeV/c 2 window around the nominal K*° mass. This includes mostly the K*° 
meson, but also an S-wave component and the interference between the S-wave and P- 
wave components. The fraction of candidates with a K*° meson present is then 



d 3 r 

Acc(fi) 



dtt 

\A S \=0 



r d 3 r 

/ Acc(fi) — dtt 
Jn ds2 



(6) 



from which the ratio f^l/f^l = 1.09±0.08 follows. Table§ summarizes all the numbers 
needed to compute the ratio of branching fractions 

B(B° ^ J/^K*°) _ 13+ o,4 ±0 50 w 

The contributions to the systematic uncertainty are also listed in Table [3] and their 
relative magnitudes are: 1.2% for the error in the efficiency ratio; 2.5% for the uncertainty 
on the transition point (a) of the Crystal Ball function; 8.6% for the parameterization of 
the upper tail of the B° peak; 3.9% for the angular correction of the efficiencies; 7.3% for 
the uncertainty on the ratio f^l / f^l an d 7.7% for the uncertainty on f d J f s . The errors 
are added in quadrature. 

Taking the value B(B° -> J/i)K* ) = (1.29 ± 0.05 ± 0.13) x 10~ 3 from Ref. |i| the 
following branching fraction is obtained, 

B{B° S -> J/^jK*°) = (4.4+°i ±0.8) x 10~ 5 . 

This value is compatible with the CDF measurement [2] and is similar to the naive quark 
spectator model prediction of Eq. (IlJ), although it is closer to the estimation in Ref. [l], 
B(B° S -> J/ipK*°) ~ 2 x B(B° d -> J/ipp°) = (4.6 ± 0.4)_x 10" 5 ._The branching fraction 
measured here is in fact the average of the B° s — > J/ipK*° and B° s — > J/ipK*° branching 
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fractions and corresponds to the time integrated quantity, while theory predictions usually 
refer to the branching fraction at t — [22]. In the case of B° s — > J/tpK*°, the two differ 
by (Aiy2r s ) 2 = (0.77 ± 0.25)%, where Ar s = T L - T H , T s = (T L + r H )/2, and T m) is 
the decay width of the light (heavy) B^-mass eigenstate. 

In conclusion, using 0.37 fb" 1 of pp collisions collected by the LHCb detector at y's = 
7 TeV, a measurement of the — > J/ipK*° branching fraction yields B(B° S — > J/ipK*°) = 
(4.4l°j ±0.8 ) x 10~ 5 . In addition, an angular analysis of the decay products is presented, 
which provides the first measurement of the K*° polarization fractions in this decay, 
giving f L = 0.50 ± 0.08 ± 0.02, /y = 0.19+^ ± 0.02, and an S-wave contribution of 
| A s | 2 = 0.07±o;o7 in a ±40 MeV/c 2 window around the K*° mass. 
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